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Abstract 


We consider an one-dimensional Schumann wave equation proposed by Marciak-Kozlowska. 
Numerical solution of that equation was obtained with the help of Mathematica. 
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Introduction 


The measured frequencies of Schuman and brainwaves are nearly the same [1]. It is worthwhile to 
point out that both calculated curves give a rather good description of the measured frequencies of 
Schuman and brain waves [2-3]. 


We consider an one-dimensional Schumann wave equation proposed by Marciak-Kozlowska. 
Numerical solution of that equation was obtained with the help of Mathematica. 


A hyperbolic equation for Schuman wave phenomena was formulated [4-5] where m is the mass 
of the neuron, fis the Planck constant, V is potential and v is the velocity propagation of the 
Schumann wave in the brain. 


Now we will obtain its numerical solution without having to recourse to Klein-Gordon equation 


as its approximation [4]. Instead, we will look for direct numerical solution and its plot using 
Mathematica 9 [6]: 


*Correspondence: Victor Christianto, Independent Researcher. Email: victorchristianto@ gmail.com 
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Concluding remarks 

A direct numerical solution of master equation corresponding to Schumann wave has been 
presented. We hope this result may be useful for further studies on the connection between 
Schuman resonance, brainwaves and other related topics. 
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Appendix: Mathematica code 

SetOptions[Plot, ImageSize->500,PlotRange->AlI,PlotPoints->nP*2,PlotS tyle- 

>{Blue, Thickness[0.01]}];{s=1/100,nP=100} 

{nN=3,l=1,11={Red,Blue,Green},12={0,1/2,1}} 

f[u_]:=2*b*a/c’2;f[u] 

eKG=D[u[x,t], {t,2} ]+a*D[u[x,t]/c, {t,1}]-D[u[x,t], {x,2}]+f[u]== 

f1C 1 [f1_]:=u[x,0]==f1 ;f1C2[f2_]:=(D[ul[x,t],t]/.t->0)==f2; 

fBC1[c_,f1_]:=(D[u[x,t],x]/.x->c)==f1; 

fBC2[d_,f2_]:=(D[u[x,t],x]/.x->d)==f2; 

{f1C 1 [f1],f1C2[f2],fBC1 [c,f1],fBC2[d,f2]}; 

paramsS={a->1,b->1,c->-1,aN->1.5};{c5=- 
5,d5=5,tF5=4,xI5=c5,xF5=d5,f15=aN*(1+Cos[2*Pi*x/d5]),f25=0,f35=0,f45=0,eKGS5=N[eKG/.params5],icS=N[{fl 
C1[f15],f1C2[f25]}/.params5],bcS=N[{fBC1[c5,f35],fBC2[d5,f45]}/.params5]} 

sol5=NDSolve[Flatten[ {eKG5,ic5,bc5 }],u,{x,xI5,xF5 }, {t,0,tF5 },MaxStepSize->s,PrecisionGoal->2] 
Do[g[i]=Plot[Evaluate[u[x,12[[i]]]/.sol5], {x,xI5,xF5 },PlotStyle- 

>{11[[]],Thickness[0.01]}], {1,1,nN}];Show[Table[g[i], (i, 1,nN}]] 

Plot3D[Evaluate[u[x,t]/.sol5], {x,xI5,xF5 }, {t,0,tF5 }, ColorFunction->Function[{ x,y}, Hue[x]],BoxRatios- 
>1,ViewPoint->{ 1,2,1},PlotRange->All,PlotPoints-> {20,20} ,ImageSize->500] 
Animate[Plot[Evaluate[u[x,t]/.sol5, {x,xI5,xF5 }],PlotRange->{-3,3}], {t,0,tF5 }, AnimationRate->0.5] 


